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Abstract.  We  report  on  recent  observations  of  AGB  stars  obtained  with  the  VLT 
Interferometer  (VLTI).  We  illustrate  in  general  the  potential  of  interferometric  measure¬ 
ments  to  study  stellar  atmospheres  and  circumstellar  envelopes,  and  demonstrate  in  particu¬ 
lar  the  advantages  of  a  coordinated  multi-wavelength  approach  including  near/mid-infrared 
as  well  as  radio  interferometry.  We  report  on  studies  of  the  atmospheric  structure  of  non- 
Mira  and  Mira  variable  giants.  We  have  used  VLTI  observations  of  the  near-  and  mid- 
infrared  stellar  sizes  and  concurrent  VLBA  observations  of  the  SiO  maser  emission.  So 
far,  this  project  includes  studies  of  the  Mira  stars  S  Ori  and  RR  Aql  as  well  as  of  the  su¬ 
pergiant  AH  Sco.  The  results  from  our  first  epochs  of  S  Ori  measurements  have  recently 
been  published  and  the  main  results  are  reviewed  here.  The  S  Ori  maser  ring  is  found  to 
lie  at  a  mean  distance  of  approximately  2  stellar  radii,  a  result  that  is  virtually  free  of  the 
usual  uncertainty  inherent  in  comparing  observations  of  variable  stars  widely  separated  in 
time  and  stellar  phase.  We  discuss  the  status  of  our  more  recent  S  Ori,  RR  Aql,  and  AH  Sco 
observations,  and  present  an  outlook  on  the  continuation  of  our  project. 
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1.  Introduction 

The  evolution  of  cool  luminous  stars,  includ¬ 
ing  Mira  variables,  is  accompanied  by  sig¬ 
nificant  mass-loss  to  the  circumstellar  envi¬ 
ronment  (CSE)  with  mass-loss  rates  of  up  to 
1(T7  -  1 0  4  Mo/year.  The  detailed  structure 
of  the  CSE,  the  detailed  physical  nature  of 
this  mass-loss  process  from  evolved  stars,  and 
especially  its  connection  with  the  pulsation 
mechanism  in  the  case  of  Mira  variable  stars, 
are  not  well  understood.  Furthermore,  one  of 
the  basic  unknowns  in  the  study  of  late-type 
stars  is  the  mechanism  by  which  usually  spher- 
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ically  symmetric  stars  on  the  asymptotic  giant 
branch  (AGB)  evolve  to  form  axisymmetric  or 
bipolar  planetary  nebulae  (PNe).  There  is  evi¬ 
dence  for  some  asymmetric  structures  already 
at  the  AGB  or  supergiant  stage  (e.g.,  Weigelt 
et  al.  1996,  1998;  Wittkowski  et  al.  1998; 
Monnier  et  al.  1999). 

Coordinated  multi-wavelength  studies 
(near-infrared,  mid-infrared,  radio,  millimeter) 
of  the  stellar  surface  (photosphere)  and  the 
CSE  at  different  distances  from  the  stellar 
photosphere  and  obtained  at  corresponding 
cycle/phase  values  of  the  stellar  variability 
curve  are  best  suited  to  improve  our  general 
understanding  of  the  atmospheric  structure, 
the  CSE,  the  mass-loss  process,  and  ultimately 
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Fig.  1.  Sketch  of  a  Mira  variable  star  and 
its  circumstellar  envelope  (CSE).  A  multi¬ 
wavelength  study  (MIDI/AMBER/VLBA)  is 
well  suited  to  probe  the  different  regions  shown 
here.  Owing  to  the  stellar  variability,  only  con¬ 
temporaneous  observations  are  meaningful. 


of  the  evolution  of  symmetric  AGB  stars 
toward  axisymmetric  or  bipolar  planetary 
nebulae.  Fig.  1  shows  a  schematic  view  of 
a  Mira  variable  star,  indicating  the  different 
regions  that  can  be  probed  by  different  tech¬ 
niques/wavelength  ranges  (VLTI/AMBER, 
VLTI/MIDI,  VLBA/maser,  ALMA). 

The  conditions  near  the  stellar  surface  can 
best  be  studied  by  means  of  optical/near- 
infrared  long-baseline  interferometry  (see, 
e.g.,  Haniff  et  al.  1995;  Wittkowski  et  al. 
2001;  Thompson  et  al.  2002;  Wittkowski  et 
al.  2004;  Woodruff  et  al.  2004;  Boboltz  & 
Wittkowski  2005;  Fedele  et  al.  2005).  The 
structure  and  physical  parameters  of  the  molec¬ 
ular  shells  located  between  the  photosphere 
and  the  dust  formation  zone,  as  well  as  of 
the  dust  shell  itself  can  be  probed  by  mid- 
infrared  interferometry  (e.g.,  Danchi  et  al. 
1994;  Ohnaka  et  al.  2005). 

Complementary  information  regarding  the 
molecular  shells  can  be  obtained  by  observ¬ 
ing  the  maser  radiation  that  some  of  these 
molecules  emit.  The  structure  and  dynamics  of 
the  CSE  of  Mira  variables  and  other  evolved 
stars  has  been  investigated  by  mapping  SiO 
maser  emission  at  typically  about  2  stellar  radii 
toward  these  stars  using  very  long  baseline 
interferometry  (VLBI)  at  radio  wavelengths 


(e.g.,  Boboltz  et  al.  1997;  Kemabll  &  Diamond 
1997;  Boboltz  &  Wittkowski  2005). 

Results  regarding  the  relationships  be¬ 
tween  the  different  regions  mentioned  above 
and  shown  in  Fig.  1  suffer  often  from  uncer¬ 
tainties  inherent  in  comparing  observations  of 
variable  stars  widely  separated  in  time  and  stel¬ 
lar  phase  (see  the  discussion  in  Boboltz  & 
Wittkowski  2005).  Both,  the  photospheric  stel¬ 
lar  size,  as  well  as  the  mean  diameter  of  the 
SiO  maser  shell  are  known  to  vary  as  a  function 
of  the  stellar  variability  phase  with  amplitudes 
of  20-50%  (Ireland  et  al.  2004;  Thompson  et 
al.  2002;  Humphreys  et  al.  2002;  Diamond  & 
Kemball  2003). 

To  overcome  these  limitations,  we  have 
established  a  program  of  coordinated  and 
concurrent  observations  at  near-infrared,  mid- 
infrared,  and  radio  wavelengths  of  evolved 
stars,  aiming  at  a  better  understanding  of  the 
structure  of  the  CSE,  of  the  mass-loss  process, 
and  of  the  triggering  and  formation  of  asym¬ 
metric  structures. 

2.  The  atmospheric  structure  of 
non-Mira  and  Mira  giants 

Fundamental  parameters,  most  importantly 
radii  and  effective  temperatures,  of  regular  cool 
giant  stars  have  frequently  been  obtained  with 
interferometric  and  other  high  angular  resolu¬ 
tion  techniques,  thanks  to  the  favorable  bright¬ 
ness  and  size  of  these  stars.  Further  param¬ 
eters  of  the  stellar  structure,  as  the  strength 
of  the  limb-darkening  effect,  can  be  stud¬ 
ied  when  more  than  one  resolution  element 
across  the  stellar  disk  is  employed.  Through 
the  direct  measurement  of  the  center-to-limb 
intensity  variation  (CLV)  across  stellar  disks 
and  their  close  environments,  interferometry 
probes  the  vertical  temperature  profile,  as 
well  as  horizontal  inhomogeneities.  Such  di¬ 
rect  limb-darkening  studies  have  been  accom¬ 
plished  for  a  relatively  small  number  of  stars 
using  different  interferometric  facilities  (in¬ 
cluding,  for  instance,  Hanbury  Brown  et  al. 
1973;  Quirrenbach  et  al.  1996;  Hajian  et  al. 
1998;  Wittkowski  et  al.  2001,  2004). 

Recent  optical  multi-wavelength  measure¬ 
ments  of  the  cool  giants  y  Sge  and  BYBoo 
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Fig.  2.  NPOI  limb-darkening  observations 
(squared  visibility  amplitude,  triple  amplitude, 
closure  phase)  of  the  M0  giant  y  Sge,  together 
with  a  comparison  to  the  best  fitting  ATLAS  9 
model  atmosphere  prediction  (squares).  For 
comparison,  the  solid  line  denotes  a  uni¬ 
form  disk  model,  and  the  dashed  line  a  fully- 
darkened  disk  model.  ATLAS  9  models  with 
variations  of  Teg  and  log  g  result  in  sig¬ 
nificantly  different  model  predictions.  From 
Wittkowski  et  al.  (2001). 


(Wittkowski  et  al.  2001)  succeeded  not  only 
in  directly  detecting  the  limb-darkening  ef¬ 
fect,  but  also  in  constraining  ATLAS  9  (Kurucz 
1993)  model  atmosphere  parameters.  Fig.  2 
shows  one  dataset  including  squared  visibil¬ 
ity  amplitudes,  triple  amplitudes,  and  closure 
phases  of  the  M0  giant  ySge  obtained  with 
NPOI,  together  with  a  comparison  to  the  best 
fitting  ATLAS  9  model  atmosphere  prediction. 
ATLAS  9  models  with  variations  of  Teg  and 
log  g  result  in  significantly  different  model 
predictions.  By  this  direct  comparison  of  the 
NPOI  data  to  the  ATLAS  9  models  alone,  the 
effective  temperature  of  ySge  is  constrained 
to  4160+100  K.  The  limb-darkening  observa¬ 
tions  are  less  sensitive  to  variations  of  the  sur¬ 
face  gravity,  and  log  g  is  constrained  to  0.9+ 1.0 
(Wittkowski  et  al.  2001).  These  results  are 
well  consistent  with  independent  estimates, 
such  as  calibrations  of  the  spectral  type. 


Fig.  3.  VLTI  limb-darkening  observations  of 
the  M4  giant  t^Phe  (Wittkowski  et  al.  2004). 


The  first  limb-darkening  observation  ob¬ 
tained  with  the  VLTI  succeeded  in  the  early 
commissioning  phase  of  the  VLTI  (Wittkowski 
et  al.  2004).  Using  the  VINCI  instrument, 
Af-band  visibilities  of  the  M4  giant  i/y  Phe 
were  measured  in  the  first  and  second  lobe 
of  the  visibility  function.  These  observations 
were  found  to  be  consistent  with  predictions 
by  PHOENIX  and  ATLAS  model  atmospheres, 
the  parameters  for  which  were  constrained 
by  comparison  to  available  spectrophotometry 
and  theoretical  stellar  evolutionary  tracks  (see 
Fig.  3).  Such  limb-darkening  observations  also 
result  in  very  precise  and  accurate  radius  esti¬ 
mates  because  of  the  precise  description  of  the 
CLV.  Future  use  of  the  spectro-interferometric 
capabilities  of  AMBER  and  MIDI  will  en¬ 
able  us  to  study  the  wavelength-dependence 
of  the  limb-darkening  effect,  which  results  in 
stronger  tests  and  constraints  of  the  model  at¬ 
mospheres  than  these  broad-band  observations 
(cf.  the  wavelength-dependent  optical  studies 
with  NPOI  as  described  above). 

For  cool  pulsating  Mira  stars,  the  CLVs  are 
expected  to  be  more  complex  than  for  non¬ 
pulsating  M  giants  due  to  the  effects  of  molec¬ 
ular  layers  close  to  the  continuum-forming 
layers.  Self-excited  hydrodynamic  model  at¬ 
mospheres  of  Mira  stars  have  been  presented 
(Hofmann  et  al.  1998;  Tej  et  al.  2003;  Ireland 
et  al.  2004,Scholz  &  Wood,  private  communi¬ 
cation).  Different  radius  definitions,  such  as  the 
Rosseland  mean  radius,  the  continuum  radius, 
or  the  radius  at  which  the  filter-averaged  inten¬ 
sity  drops  by  50%,  may  result  in  different  val- 
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ues  for  the  same  CLV.  On  these  topics,  see  also 
Scholz  (2003).  However,  interferometric  mea¬ 
surements  covering  a  range  of  spatial  frequen¬ 
cies  can  directly  be  compared  to  CLV  predic¬ 
tions  by  these  model  atmospheres  without  the 
need  of  a  particular  radius  definition.  At  pre¬ 
maximum  stellar  phases,  when  the  temperature 
is  highest,  the  broad-band  CLVs  are  less  con¬ 
taminated  by  molecular  layers,  and  different 
radius  definitions  agree  relatively  well  (Scholz 
&  Wood,  private  communication). 

A''- band  VINCI  observations  of  the  proto¬ 
type  Mira  stars  o  Cet  and  R  Leo  have  been  pre¬ 
sented  by  Woodruff  et  al.  (2004)  and  Fedele  et 
al.  (2005),  respectively.  These  measurements 
at  post-maximum  stellar  phases  indicate  in¬ 
deed  A"- band  CLVs  which  are  clearly  different 
from  a  uniform  disk  profile  already  in  the  first 
lobe  of  the  visibility  function.  The  measured 
visibility  values  were  found  to  be  consistent 
with  predictions  by  the  self-excited  dynamic 
Mira  model  atmospheres  described  above  that 
include  molecular  shells  close  to  continuum¬ 
forming  layers. 

3.  Joint  VLTI/VLBA  observations  of 
the  Mira  star  SOri 

We  started  our  project  of  joint  VLTI/VLBA 
observations  of  Mira  stars  in  December 
2002/January  2003  with  coordinated  near- 
infrared  K -band  VLTI/VINCI  observations  of 
the  stellar  diameter  of  the  Mira  variable  S  Ori 
and  quasi-simultaneous  VLBA  observations  of 
the  43.1  GHz  and  42.8  GHz  SiO  maser  emis¬ 
sions  toward  this  star  (Boboltz  &  Wittkowski 
2005).  We  obtained  in  December  2004/January 
2005  further  concurrent  observations  including 
mid-infrared  VLTI/MIDI  observations  to  probe 
the  molecular  layers  and  the  dust  shell  of  S  Ori, 
and  new  epochs  of  VLBA  observations  of  the 
43.1  GHz  and  42.8  GHz  SiO  maser  rings. 

The  December  2002/January  2003  repre¬ 
sent  the  first-ever  coordinated  observations  be¬ 
tween  the  VLTI  and  VLBA  facilities,  and  the 
results  from  these  observations  were  recently 
published  in  Boboltz  &  Wittkowski  (2005). 
Analysis  of  the  SiO  maser  data  recorded  at  a 
visual  variability  phase  0.73  show  the  average 
distance  of  the  masers  from  the  center  of  the 


Lightcurve  of  S  Ori 


Julian  Date  -  2450000 


Fig.  4.  Lightcurve  of  S  Ori  together  with  the 
epochs  of  our  joint  VLTI/VLBA  measurements 
obtained  so  far.  Note  that  the  y-axis  is  given 
with  increasing  V  magnitude,  i.e.  the  stellar 
maximum  is  at  the  bottom  and  stellar  min¬ 
imum  at  the  top.  The  study  of  S  Ori  was 
started  in  ESO  period  P70  (Dec.  2002/Jan. 
2003)  including  near-infrared  K- -band  VINCI 
and  VLBA/SiO  maser  observations  (Boboltz  & 
Wittkowski  2005).  In  December  2004/January 
2005,  we  obtained  concurrent  mid-infrared 
VLTI/MIDI  and  VLBA/SiO  maser  observa¬ 
tions. 


distribution  to  be  9.4  mas  for  the  v  =  1,/  = 
1-0  (43.1  GHz)  masers  and  8.8  mas  for  the 
v  =  2,  /  =  1  -  0  (42.8  GHz)  masers.  The 
velocity  structure  of  the  SiO  masers  appears 
to  be  random  with  no  significant  indication 
of  global  expansion/infall  or  rotation.  The  de¬ 
termined  near-infrared,  K -band,  uniform  disk 
(UD)  diameters  decreased  from  ~  10.5  mas  at 
phase  0.80  to  ~10.2  mas  at  phase  0.95.  For  the 
epoch  of  our  VLBA  measurements,  an  extrap¬ 
olated  UD  diameter  of  @fjD  =  10.8  +  0.3  mas 
was  obtained,  corresponding  to  a  linear  radius 
of  R*d  =  2.3  ±  0.5  AU  or  R*p  =  490  ± 
115  Rq.  The  model  predicted  difference  be¬ 
tween  the  continuum  and  K -band  UD  diame¬ 
ters  is  relatively  low  in  the  pre-maximum  re¬ 
gion  of  the  visual  variability  curve  as  in  the 
case  of  our  observations  (see  above).  At  this 
phase  of  0.73,  the  continuum  diameter  may  be 
smaller  than  the  K -band  UD  diameter  by  about 
15%  (Ireland  et  al.  2004).  With  this  assump¬ 
tion,  the  continuum  photospheric  diameter  for 
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Fig.  5.  First-ever  coordinated  observations  be¬ 
tween  ESO’s  VLTI  and  NRAO’s  VLBA  facili¬ 
ties:  43.1  GHz  SiO  maser  emission  toward  the 
Mira  variable  S  Ori  measured  with  the  VLBA, 
together  with  the  near-infrared  diameter  mea¬ 
sured  quasi  simultaneously  with  the  VLTI  (red 
stellar  disk).  From  Boboltz  &  Wittkowski 
(2005). 


the  epoch  of  our  VLBA  observation  would  be 
©Phot(VLBA  epoch,  phase  =  0.73)  »  9.2  mas. 
Our  coordinated  VLBA/VLTI  measurements 
show  that  the  masers  lie  relatively  close  to  the 
stellar  photosphere  at  a  distance  of  ~  2  photo- 
spheric  radii,  consistent  with  model  estimates 
by  Humphreys  et  al.  (2002)  and  observations 
of  other  Mira  stars  by  Cotton  et  al.  (2004). 
The  new  2004/2005  VLTI  and  VLBA  data  are 
currently  being  reduced  and  analyzed. 


4.  Outlook 

We  concentrate  on  a  few  stars  in  order  to  un¬ 
derstand  the  CSE  for  a  few  sources  in  depth. 
In  addition  to  the  S  Ori  data  described  above, 
we  have  to  date  VLTI/MIDI  observations  of 
the  supergiant  AH  Sco,  and  of  the  Mira  star 


RR  Aql,  as  well  as  concurrent  VLBA  obser¬ 
vations  for  each  of  these  targets/epochs.  These 
data  are  currently  being  analyzed. 

Further  studies  will  aim  at  including  more 
detailed  near-infrared  studies  of  the  stellar 
atmospheric  structure  (close  to  the  photo¬ 
sphere)  employing  VLTI/ AMBER,  concurrent 
with  VLTI/MIDI  and  VLBA  observations  as 
discussed  above.  Making  use  of  the  spectro- 
interferometric  capabilities  of  AMBER,  and 
of  the  closure-phase  information,  these  studies 
can  in  principle  reveal  horizontal  surface  inho¬ 
mogeneities  (see  Wittkowski  et  al.  2002). 
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